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A more precise determination of the effective fine structure constant a e s(E) is mandatory for confronting data 
from future precision experiments with precise SM predictions. Higher precision would help a lot in monitoring 
new physics by increasing the significance of any deviation from theory. At a future e + e~-collider like the ILC, 
as at LEP before, a e g(E) plays the role the static zero momentum a = a e fr(0) plays in low energy physics. 
However, by going to the effective version of a one loses about a factor 2 x 10 2 at E — to 10 5 at E = Mz 
in precision, such that for physics at the gauge boson mass scale and beyond a c s(E) is the least known basic 
parameter, about a factor 20 less precise than the neutral gauge boson mass Mz and by about a factor 60 less 
precise than the Fermi constant Gf- Examples of precision limitations are a e ff(ro M ) which limits the theoretical 
precision of the muon anomalous magnetic moment a M and a e ff(Mz) which limits the accuracy of the prediction 
of the weak mixing parameter sin 2 0/ and indirectly the upper bound on the Higgs mass mu- An optimal 
exploitation of a future linear collider for precision physics requires an improvement of the precision of a e g(E) 
by something like a factor ten. We discuss a strategy which should be able to reach this goal by appropriate 
efforts in performing dedicated measurements of (Jhadronic in a wide energy range as well as efforts in theory and 
in particular improving the precision of the QCD parameters a s , m c and mj by lattice QCD and/or more precise 
determinations of them by experiments and perturbative QCD efforts. Projects at VEPP-2000 (Novosibirsk) and 
DANAE/KLOE-2 (Frascati) are particularly important for improving on a e g(Mz) as well as a or r(nv). Using 
the Adler function as a monitor, one observes that we may obtain the hadronic shift Aa had (Af§) as a sum 
Aa h 5 a ) d (-s ) data + Aa^ d (s , M|) pQCD where the first term includes the full non-perturbative part with the 
choice so = (2.5 GeV) 2 or larger. In such a determination low-energy machines play a particularly important role 
in the improvement program. We present an up-to-date analysis including the recent data from KLOE, SND, 
CMD-2 and BABAR. The analysis based on e + e"-data yields Aa^(M|) = 0.027593 ± 0.000169 [a _1 (M|) = 
128.938 ±0.023] (splitting with s = (10 GeV) 2 to reduce dependence on m c ), Aa^ d (M§) = 0.027607 ± 0.000225 
[« _1 (Mi) = 128.947± 0.035] (standard approach) , and a£ ad = (692.1 ±5.6) xl0~ 10 . The continuation of a e fi (E) 
from the Z mass scale to ILC energies may be obtained by means of perturbative QCD. We emphasize the very 
high improvement potential of the VEPP-2000 and DANAE/KLOE-2 projects. 



1. INTRODUCTION 

Of all non-perturbative hadronic effects enter- 
ing the predictions of a host of electroweak pre- 
cision observables, the main effect enters via the 
effective fine-structure "constant" a(E) |1I2| . It 
is the non-perturbative hadronic contribution to 
charge screening by vacuum polarization (VP) 
which limits its precision. The value of a(E) 
is of interest for a wide range of energies, well 
known examples are a(Mz) or a(m^) where the 
latter accounts for the leading hadronic contri- 
bution in the muon anomaly = (g — 2) M /2 



(see |3I4I5I6I7I8I9I1U) for some recent e + e~ data 
based evaluations). While electroweak effects 
like lepton contributions are calculable in per- 
turbation theory, for the strong interaction part 
(hadron- and/or quark- contributions) perturba- 
tion theory fails and a dispersion integral over 
e + e~-data encoded in 

o-(e~*"e — — > 7* — > hadrons) 

= i + _ ; + \ 

cr(e+e — > 7* — > ) 

provides a reliable approach to estimate the non- 
perturbative effects. Errors of data directly en- 
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s' = M} (1 - k) [k = E,/^] 



Energy scan Photon tagging 

Figure 1. The two schemes of measuring 
cr(e + e _ — > hadrons) . 

ter the theoretical uncertainties of any prediction 
depending on a(-E) at non-zero E. Evaluations 
of the above-mentioned dispersion relation have 
developed into an art of getting precise results 
still often from old cross-section measurements 
of poor precision. What is needed is a reduc- 
tion of the present error by an order of magnitude 
within the next ten years. This is a new challenge 
for precision experiments on cr(e + e~ — > hadrons) 
such as ongoing experiments KLOE, BABAR, 
and Belle which measure <7h a dronic via radiative 
return or photon tagging (see Fig^l |11I12I13| . 
For future precision experiments, in many cases, 
we will need to know the running a om very pre- 
cisely, desirably at the per mill level. Note that 
corrections are large and steeply increasing at low 
E, as may be learned from Fig [5] 

An immediate question might be: why not 
measure a e s(E) directly, like the QCD run- 
ning coupling a s (s)? The problem is that any 
measurement requires a normalizing process like 
Bhabha (Fig|3Jl which itself depends on a c g(t) 
and a e g(s). In fact one is always measuring some- 
thing like 

r(E) oc (a cS (s)/a cS (t)) 2 , (1) 

with t = — i (s — Am 2 e ) (1 — cos#). Unless one 
is able to measure essentially at zero momentum 
transfer (~ angle), because of the steep rise of 
a e fi(E) at low energies a large fraction of the ef- 
fect which we would like to determine drops out, 
especially the strongly rising low-energy piece, 
which includes substantial non-perturbative ef- 
fects (see also HU). 

It should be noted that not only a e g(E) at 
high energies is of interest. As a logarithmi- 
cally increasing function the problems show up 
at relatively low scales like for a c g(m^), which 
determines the leading uncertainty of the muon 
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Figure 2. The running of a(E) = a(0)/(l — Aa) 
in terms of Aa. The "negative" E axis is chosen 
to indicate space-like momentum transfer. The 
vertical bars at selected points indicate the uncer- 
tainty. In the time-like region the resonances lead 
to pronounced variations of the effective charge 
(shown in the p — oj and <f> region). 

anomaly (see Fig. 0J. Another not so well 
known example is a c g(M pmton ) at the proton 
mass scale which enters /3-decay and affects the 
determination of the CKM element V u d |15I16| . 

In spite of the fact that the discussion presented 
here is not really new, the substantial amount of 
new data in the low energy region from KLOE, 
SND and CMD-2 and from BABAR leads to a no- 
ticeable change of the error profile of a e g(E) at 
different scales and it provides a strong motiva- 
tion for the upcoming VEPP-2000 Q7| and future 
DANAE/KLOE-2 [H] projects to continue pre- 
cise measurements of Uhadronic m the region from 
threshold up to 2.5 GeV. 

2. a(M z ) IN PRECISION PHYSICS 

At higher energies for all processes which are 
not dominated by a single one photon exchange, 
a e s(E) enters in a complicated way in observab- 
les and cannot be measured in any direct way. In 




Figure 3. Running a in Bhabha scattering as a 
normalization process. The "blobs" are VP inser- 
tions into the tree level diagrams. 



Myy and to the weak mixing parameter sin 2 9 * 



as 



SMw 
M w 
S sin 2 6 / 
sin 2 f 



sin Gv 



2 cos 2 Gw — sin 2 Qi 
cos 2 Q f 



SAa ~ 0.23 8Aa 



cos 2 0/ — sin 2 f 



SAa ~ 1.54 SAa 



and to all kinds of observables which depend 
on these parameters. In particular the indirect 
bounds on the Higgs mass obtained from elec- 
troweak precision measurements are weakened 
thereby. 




Figure 4. Leading hadronic contribution to g M — 2 



3. EVALUATION OF a{M z ) 

The non-perturbative hadronic shift Aa n ^ d (s), 
due to the 5 light quark flavors, can be evaluated 
in terms of a(e + e~ — ► hadrons) data via the well 
known dispersion integral (see [3] and references 
therein): 



the SM many parameters are interrelated as im- 
plied by the mass generating Higgs mechanism. 
A nice example are parameters showing up in 
typical four fermion and vector boson processes. 
Unlike in QED and QCD the SM is a sponta- 
neously broken gauge theory which leads to a new 
kind of parameter interdependencies, which have 
been tested at LEP with high accuracy. Besides 
the fermion masses and mixing parameters the 
SM has only 3 independent parameters the gauge 
couplings: g and g' and the Higgs vacuum ex- 
pectation value v, which may be fixed from the 
most precisely measured quantities, namely, a, 
Gfj, and Mz- All other dependent parameters are 
then predictions which can be tested and provide 
a monitor for new physics. 

The impact of the hadronic uncertainties of 
a e g(Mz) in physics of the heavy gauge bosons 
M\y and Mz are well known and for a more de- 
tailed discussion I refer to my earlier articles pQ. 
As mentioned before, in place of a, which is 
known to a precision — ~ 3.6 x 10~ 9 , the ef- 
fective a c g(Mz) is needed and non-perturbative 
hadronic effects reduce its precision to S "^^-j ~ 

1.6 -7- 6.8 x 10~ 4 , while the other two basic pa- 
rameters are known with much better accuracy 

^ - 8.6 x 1CT 6 and ~ 2.4 x 1(T 5 . The un- 

Ct m Mz 

certainty in a c g(Mz) carries over to the W mass 
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where 



r(°)(e+e- 



hadrons) 



(2) 



(3) 



3s 



A compilation of the data is shown in Fig. [S] 
For the evaluation at Mz — 91.1876 GeV we take 

• R(s) data up to y/s = E cut = 5.2 GeV and for 
T resonances region between 9.6 and 13 GeV, 

• perturbative QCD pErT] from 5.2 to 9.6 GeV 



and for the high energy tail above 13 GeV, as 
recommended in |21j . As a result we obtain 



hadrons 



(Mf) = 0.027607 ± 0.000225 



128.947 ± 0.035. 



(4) 



Contributions from various energy regions and 
the origin of the errors are shown in Fig. HJ 

As we will explain below Aa^ d (M z ) , or more 
generally a e g(E) for E > 2.5 GeV, may be eval- 
uated in a different way, by exploiting pQCD as 
much as possible in a well controlled fashion via 
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Figure 5. A compilation of the presently available 
experimental hadronic e + e~ -annihilation data 

the Adler function. In this Adler function ap- 
proach the complete non-perturbative part may 
be evaluated at y/so — 2.5 GeV where 

Aa^j(-so) = 0.007364 ±0.000101. 

and Aa^j(Mf) — Aa^ d (-s ) is reliably calcu- 
lable using pQCD. The profile of Aa llad (—so) is 
shown in Fig. for more details we refer to Table 
below. 

4. A LOOK AT THE e+e' DATA 

In order to learn where results have to be im- 
proved we briefly have a closer look at the ex- 
isting e + e~ — > hadrons cross-section data. Since 
our analysis in 1995 data from MD1 [22], BES- 
II HS| and from CMD-2 |23 have lead to a sub- 
stantial reduction in the hadronic uncertainties 
on Aa had and a^ ad . More recently KLOE (25] . 





Figure 6. Aaj^ d (M|): contributions (left) and 
errors 2 (right) from different regions 




Figure 7. Aal^L(— So): contributions (left) and 
errors 2 (right) from different regions 

SND I2H1 and CMD-2 [27] published new mea- 
surements in the region below 1.4 GeV. Unfor- 
tunately the agreement between the different ex- 
periments is not very satisfactory. Nevertheless, 
the progress is substantial, the low energy domain 
is no longer essentially dominated by one exper- 
iment (CMD-2 2003) and the errors of the com- 
bined data are reduced noticeably and existing 
problems can be settled by ongoing experiments. 
In the past few years also more data on the purely 
neutral channels 7r°7, 7r°7r°7, 777, 7r°?77, 7/7, ujtt° 
mainly from SND [2HI and CMD-2 |2S] have been 
included. Substantial improvement on the u> and 
resonances, on Air and 5tt channels and a study 
of r77r + 7r~, ujt: + tt~ |30| contributed to the im- 
provement. An important conclusion from the 
CMD-2 study |2H] is that the u>- and ^-resonance 
decays completely saturate the cross sections of 
these channels. Consequently, there are no unac- 
counted contributions to the R value from such 
neutral channels (see [2] for a more detailed ex- 
pert summary). Further data included 3tt, 4ir, 
KK and pp channels and come from SND, CMD- 
2, BES-II and CLEO |32I33I34I35I5SI57| 

As in other more recent evaluations, modes 
not measured directly have been included by es- 
timating them using isospin relations and known 
branching fractions of decay modes. Included are 



(see also 0E]) (1) 77+77-371-° ~ [2(77+77-) tt° - 
?77r + 7r")]/2 + 777r + 7r~ -BR(r) -> 2ir°), (2) 777777 not 
included already in the 5ir modes, (3) uj(ui — > 
7r°7)7r7r including the aj7r°7r° mode via isospin 
(= (ujir+ir- + W7r°7r°[= l/2w7r+7r-]) ■ BR(w -> 
7r°7), (4) unseen KKirir modes using e + e~ — > 
-ft^A" data of DM1, which amounts to taking 2 
times the total for the modes i^°A' (7r7r) and 
K+K-{niif (=2(K° S X - K° S K° L - K^K^ir* - 
K+R-tt )), K+K-jt+tt- measured by DM1 and 
BABAR is evaluated from the data and hence is 
to be subtracted from the total KKirir contribu- 
tion, (5) missing isospin modes in 6ir: 77 + 77 _ 47r° 
(~ 0.093 • 3(77+77-) + 0.031 • 2(77+77-) 2tt°), (6) 
P-f (P — 7r°,ry) (note 7r°7r°7 is a version of 
urn already included), (7) K + K~ir° (was miss- 
ing in non-DM2 R compilations) also accounting 
for KlK%Tr°[= K+R-ir }. Channel measured 
for the first time include K + K~ 2(77+77-) and 
2(K+K-) (BABAR). For a useful review on the 
e+e~ data we refer to [2U ( U P to 2003). 

At higher energies data are particularly prob- 
lematic in the region between 1.4 and 2.5 GeV. 
Fortunately, a new set of measurements is avail- 
able now from radiative return experiments at 
BABAR for the exclusive channels e+e~ — » 

TT+TT-TT , TT+TT-TT+TT-, K+R-JT+T:- , 2(K+K~), 

3 (vr+O, 2(77+77-77°) and K+K^ic+Tr ). 
These data cover a much broader energy interval 
and extend to much higher energies than previous 
experiments. The compilation of all data in this 
range is shown in Fig. [5] Some of the early exper- 
iments measured exclusive processes channel by 
channel, up to 2.5 GeV including about 23 chan- 
nels (see FiglHJ and others performed inclusive 
measurements for R(n > 2) to which the two- 
body channels have to be added. The latter drops 
to a small contribution above 1.4 GeV. This is in 
favor of an inclusive strategy and helps to sepa- 
rate leptonic two prong events from the hadronic 
channels dominated by n > 2 events. In view 
of the many channels an inclusive measurement 
seems to be more feasible at a precision of about 
1% which should be attempted in the improve- 
ment program. Table Ogives a more detailed pic- 
ture about the relevance of the various modes for 
the contribution to the range 2Mr- < E < 2 GeV. 
In our analysis we are working throughout with 



Table 1 

Contributions to ajj ad and Aa|^j(- sq) from the 
energy region 2M K < E < 2 GeV. X* = 
7r °7) 5 i SO =evaluated using iospin relations. 



channel X j a* % | Aa x % 



7T°7 


0.04 


0.04 


0.00 


0.03 




11.99 


11.66 


1.59 


9.64 




9.22 


8.98 


1.25 


7.55 


VI 


0.45 


0.44 


0.05 


0.30 




19.27 


18.75 


3.79 


22.93 


27T + 27T~ 


13.99 


13.62 


2.80 


16.92 


7T + 7T~37r iBO 


1.17 


1.14 


0.26 


1.56 


27r+27r~7T° 


1.94 


1.88 


0.43 


2.60 


7r+7r _ 47r° iao 


0.08 


0.08 


0.02 


0.12 




0.26 


0.25 


0.05 


0.31 


27r+2ir _ 2ir° 


1.70 


1.65 


0.42 


2.54 


37r+37r~ 


0.32 


0.31 


0.08 


0.49 




0.77 


0.75 


0.13 


0.78 


K + K~ 


21.99 


21.39 


2.64 


15.94 


K%Kl 


13.17 


12.82 


1.49 


8.99 




0.09 


0.08 


0.02 


0.12 


K+K-tt 


0.35 


0.34 


0.08 


0.49 


k° s kU° iso 


0.35 


0.34 


0.08 


0.49 




1.08 


1.05 


0.25 


1.49 




1.08 


1.05 


0.25 


1.49 


K + K~ 7T +77 ~ 


1.08 


1.05 


0.28 


1.70 


KKlTTT iBO 


2.22 


2.16 


0.54 


3.23 


pp 


0.07 


0.07 


0.02 


0.12 


nn 


0.08 


0.07 


0.02 


0.13 


— > missing 


0.03 


0.03 


0.00 


0.02 


sum 


102.78 


100.00 


16.54 


100.00 


tot [sum in %] 


692.00 


[14.85] 


73.65 


[22.46] 



renormalized bare cross sections © ■ For older 
data the missing renormalizations are performed 
accordingly. Not performing these a posteriori ra- 
diative corrections would lead to a value for a'^ ad 
which is about 1 a higher. Final state photon 
radiation for the low energy dominating 7r + 7r~ 
channel is included as given by scalar QED. In 
fact this procedure is less model dependent than 
it looks like in first place. The reason is that ex- 
periments have subtracted radiative events using 
the same model, while in fact at least the virtual 
hard photon effects are included in the measured 
cross sections, as virtual effects cannot be elimi- 
nated by cuts |41| . 

An intresting reconsideration of data in the 
J/ip resonance region |42| demonstrates very 
good agreement between BES and much older 
Crystal Ball 02! results and in fact taking into ac- 
count only the two data sets leads to a reduction 
of the uncertainty in this region. 



Table 2 

Contributions for Aaj^(M§) x 10 4 (direct integration method) and Aa[ 5 ^(- s ) x 10 4 (non-perturbative 
part in the Adler function method), with relative (rel) and absolute (abs) error in percent. 



Energy range 




(Aff) x 10 4 


rel [%} 


abs [%] 


A°C(-*o) x 104 


rel [%} 


abs [%] 


p,w (E < 2M K ) 


36.23 


13.1] (0.24) 


0.7 


1.1 


33.29 [ 45.2] (0.22) 


0.7% 


4.6 % 


2M K < E < 2 GeV 


21.80 


7.9] (1.33) 


6.1 


34.9 


16.44 [ 22.3] (0.83) 


5.0 % 


67.7 % 


2 GeV < E < M J/4 , 


15.73 


5.7](0.88) 


5.6 


15.4 


7.91 [ 10.7] (0.44) 


5.6 % 


19.3 % 


% < E < My 


66.95 


24.3] (0.95) 


1.4 


18.0 


13.94 [ 18.9] (0.29) 


2.0 % 


8.1 % 


My < E < E cut 


19.69 


7.1](1.24) 


6.3 


30.4 


0.96 [ 1.3](0.06) 


6.2 % 


0.4 % 


E cut < E pQCD 


115.66 


41.9](0.11) 


0.1 


0.3 


1.09 [ 1.5](0.00) 


0.1 % 


0.0 % 


E < E cut data 


160.41 


58.1](2.24) 


1.4 


99.7 


72.55 [ 98.5](1.01) 


1.4% 


100.0 % 


total 


276.07 


100.0] (2.25) 


0.8 


100.0 


73.64 [100.0](1.01) 


1.4% 


100.0 % 
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Figure 8. Thresholds for exclusive multi particle 
channels below 2 GeV 



Substantial progress is expected in this region 
by the future VEPP-2000 [TZ\ at Novosibirsk and 
a possible new facility at Frascati JSJ- Not only 
the effective fine structure constant can be sub- 
stantially improved, also the hadronic contribu- 
tion to the muon g — 2 may be improved (see 
TabEJ 

The present error profile of the various quanti- 
ties of interest is compared in FigllOl 



5. Aa had AND THE ADLER FUNCTION 

The Adler function is an ideal tool for disentan- 
gling perturbative from non-perturbative effects 
in the Euclidean region. It is defined by 



n -s = — a —Aa (s = - (12tt 2 ) s — 

a ds ' ds 



and thus for the hadronic part we may write 

oc 

Ry(s) 



ds ■ 



4ffi2 



(s + Q 2 



(5) 



While the time-like function R 7 (s) is calculable 
in pQCD only by refering to quark-hadron dual- 
ity, D 1 {Q 2 ) is a smooth simple function both in 
terms of hadrons (dispersion integral over physi- 
cal cross sections) and in terms of quarks (pQCD) 
such that the validity of pQCD can be exam- 
ined directly. A comparison of the experimental 
vs. the pQCD Adler function Fig. II II shows that 
pQCD in the Euclidean region works very well for 
yQ2>2.5 GeV glj. The point is that at energies 
above that "threshold" massive QCD works in 
any renormalization scheme provided mass effects 
are taken into account correctly 1 . The breakdown 
of pQCD is obviously due the fact that we are us- 
ing the MS running coupling a^ s (s) which grows 



1 Note that the curvature of D 1 {Q 2 ) is almost solely due 
to the quark masses. In massless QCD in the MS scheme 
D 7 (Q 2 ) is essentially a constant depending on the number 
of flavors Nf. 



Table 3 

Contributions to ajj ad x 10 10 with relative (rel) and absolute (abs) error in percent. 



Energy range 


< d [%]( 


error) x 10 10 


rel [%} 


abs [%] 


p,u (E< 2M K ) 


538.33 


[ 77.8] (3.65) 


0.7 % 


42.0 % 


2M K < E < 2 GeV 


102.31 


[ 14.8](4.07) 


4.0 % 


52.1 % 


2 GeV < E < 


22.13 


[ 3.2] (1.23) 


5.6 % 


4.8 % 


M JM < E < M T 


26.40 


[ 3.8] (0.59) 


2.2 % 


1.1 % 


M r < E < E cut 


1.40 


[ 0.2] (0.09) 


6.2 % 


0.0 % 


E cut < E pQCD 


1.53 


[ 0.2](0.00) 


0.1 % 


0.0 % 


E < E cut data 


690.57 


[ 99.8] (5.64) 


0.8 % 


100.0 % 


total 


692.10 


100.0](5.64) 


0.8 % 


100.0 % 



to oo at y/s = Aq^ d 2 . In any case it looks con- 
vincing to calculate 



AaS d (M 



A«£!d(-*o) 



(lain 



(6) 



A«E(-^I)-A< a J d (- S0 ) 



(5) 



A4 5 2i(^I)-Aa^ d (-Ml) 



(5) 



pQCD 



pQCD 



where the first term only has to be calculated 
from the experimental data. The large pertur- 
batiye second term now is rather sensitive to a 
precise knowledge of the QCD parameters a s , m c 
and mj, however. For our pQCD evaluation we 
take the QCD parameters: 

a s (M z ) = 0.1189 ± 0.0010 @7], 
and quark masses (in GeV): 



masses 


non-lattice 


lattice only 


m c (m c ) 
fh b (fh b ) 


1.24 ± 0.09 M8J 
4.20 ± 0.07 48 


1.30 ± 0.03 ± 0.20 K9J 
4.20 ± 0.10 ± 0.10 51 



as reviewed in [501 . Lattice results now provide 
important cross checks of the phenomenological 

2 One could attempt to go to an infrared finite renormal- 
ization scheme e.g. by redefining in pQCD 

aF(s) 



-scheme 



(-) 



i - 



a 



MS 



{s)/c 



where c a constant O(l). Now a c s scheme ^ s ^ g Des to the 
finite value c (of your chioce) as yfs — > Aq^ D , which ap- 
parently would extend the validity of the pQCD result to 
lower energies. But this would be at the expense of a huge 
scheme dependence. The point is that a s is not by itself an 
observable and unless one precisely specifies an observable 
which defines it, its meaning remains unclear. Also, this 
kind of regularization "by hand" introduces an infrared 
fixed point in the /3-function which could be in conflict 
with confinement expectations. For a different point of 
view see 1451 and references therein. 



sum rule results which we will actually use. 

In this approach (see pQ for a detailed discus- 
sion of the parameter dependences and error es- 
timates) we obtain 



Aa£i(-*o) 



data 



0.007364 ±0.000101 (7) 



and together with the second pQCD term we ar- 
rive at 



Aa n a } d (-^i) 



0.027478 ±0.000144 
128.954 ± 0.020. 



a\-M\) 
Finally, for the third term of I© we obtain 

A = 0.000038 ±0.000002, 
such that our final result reads 



Aa£,(M!) 



" l (M 



0.027547 ±0.000144 
128.949 ± 0.020. 



(8) 



We notice that this result already has a substan- 
tially lower error than the value obtained via the 
direct dispersion integral. Errors of the perturba- 
tive part have been taken to be 100% correlated 
(worst case). The sensitivity to the charm mass 
in particular is substantial: 



parameter 


range 


pQCD uncertainty 


ct s 


0.1179 ... 0.1199 


0.000017 


m c 


1.15 ... 1.33 


0.000078 


m b 


4.13 ... 4.27 


0.000008 



but it can be reduced by choosing at higher scale 
like y/so — 10 GeV to diminish the m c depen- 
dence. The value obtained is given in the abstract 
as my best conservative estimate. 

The conclusion of our analysis is that while 
in the time-like approach pQCD works well only 
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Figure 9. Status of exclusive and inclusive mea- 
surements in the most problematic region 

in the "perturbative windows" 3.00 — 3.73 GeV, 
5.00-10.52 GeV and 11.50- oo, in the space-like 
approach the plot of the Adler function shows 
that pQCD works well for all Q 2 = -q 2 > 2.5 
GeV. 

For the future, in particular in view of ILC re- 
quirements, one should attemt to improve the de- 
termination of the effective a cm by a factor 10 in 
accuracy. What may we expect to be realistic? 
We assume that dedicated cross section measure- 
ments are possible at the 1% level in the relevent 
energy regions. One then obtains the following 
estimates: 

• with the direct integration of the data, and 
cuts applied as above, 58% of the hadronic contri- 
bution to a(Mz) is obtained from data and 42% 
from pQCD. Given Aa| i 5 a ) d data x 10 4 = 160.41 ± 
2.24 (1.4% error) assuming a 1% overall accuracy 
the error would be ±1.60. However, assuming 
that different independent experiments are per- 
formed at 1% accuracy for each region (divided 
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Figure 10. Comparison of error profiles between 
M 5 2l(^f), A4 5 a } d (- So ) and a. 



up as in Tab. |2J and adding up errors in quadra- 
ture, one gets ±0.83. The improvement factor 
from the data thus would be 2.7. The theory part 

Aa had PQCD x 104 = H5.66±0.11 already now has 
an 0.1% accuracy and would not be required to 
be improved. 

• Using integration via the Adler function we 
have a 26% contribution from data and a 74% 



(5) data 
had 



pQCD. The experimental part is Aa 
10 4 = 73.64 ± 1.01 (1.5% accuracy) and a 1% 
overall accuracy would reduce the error to ±0.74. 
Again, assuming that independent 1% accuracy 
measurements are possible for each region (di- 
vided up as in Tab. |5J) and combining errors in 
quadrature would yield ±0.40. Again the im- 
provement factor is essentially as above 2.8 but 
5.6 with respect to the present direct method 
which is usually adopted. The important point, 
however, is the very different error profile Fig. 1101 
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Figure 11. "Experimental" Adler- function versus 
theory (pQCD + NP) in the low energy region (as 
discussed in |44|). Note that the error includes 
both statistical and systematic ones, in contrast 
to Fig. where only statistcal errors are shown. 

in the two approaches. In the Adler function ap- 
proach errors may be reduced to a large extent 
alone by low energy machines below 2.5 GeV. 

• One important draw back is that the large 
pQCD part Aa 1 ( 1 5 a ) d pQCD x 10 4 = 201.83 ± 1.03 
uses pQCD down to much lower energies, and 
the parameter uncertainties become much more 
severe leading to a 0.5% accuracy "only". Here 
an improvement by a factor 5 would be desirable. 
There has been steady progress in the past and 
we have no doubt that much improvement will be 
possible in coming years. 

6. CONCLUSION 

The analysis presented above suggests that 
an improvement by a factor 5 in the error of 
a e g(Mz) may be realistic. The strategy for 
reaching this goal: 

• Adopt the Adler function approach to mon- 
itor and control non-perturbative strong interac- 
tion effects, i.e., write Aa^ d (M^) in the form 

and determine Aa^^- sq) from data using the 
dispersion relation J5J. 

• The determination of Actwj( — so) requires a 
dedicated program of a(e + e~ — > hadrons) cross 
section measurements especially at low energy 
machines attempting a precision of 1% in the to- 
tal cross section. As much as possible the inclu- 



sive method should be pushed as it seems to be 
simpler and easier to reach the 1% accuracy. De- 
tailed Monte Carlo simulation and detector stud- 
ies are necessary to clarify the possibilities. 

• Such experiments have to be accompanied 
by QED and SM calculations of processes like 
Bhabha scattering, /i-pair production and 7r-pair 
production at least at the level of complete two- 
loop QED plus resummations of leading higher 
order effects, and leading weak effects. 

• The evaluation of the missing pieces needed 
to obtain e.g. Aa|^ d (Af|) may be performed us- 
ing pQCD. This 74% piece must be evaluated 
at a precision at least at the 0.5% level. The 
four-loop pQCD calculation of the Adler function 
should be extended to include mass effects. Much 
more important is a substantial improvement in 
the precise determination of the QCD parameters 
a s , m c , and, to a lesser extent, mj. Here lattice 
QCD must play a crucial role (see e.g. and 
references therein). 

• Obviously, to reach this goal requires a big 
effort especially on the experimental side. Of 
course the simple straightforward direct integra- 
tion approach further will be applied and will pro- 
vide an important cross check for the evaluation 
based on the splitting (JSJ. However, to reach the 
same precision using this standard method would 
require much more experimental effort also at 
higher energies as may be concluded from Fig. llUI 

In the near future progress will be possible 
by the radiative return experiments at KLOE, 
BABAR and Belle g0] up to ^/i ~ 3 GeV and 
by a new energy scan experiment at VEPP-2000 
up to y/s — 2 GeV. Further improvements of 
the accuracy of the R measurements in the range 
3 GeV < y/s < 5 GeV is also expected from 
CLEO-Cjnnj and from the r-charm factory with 
BES-III jH] in Beijing. At Frascati a new facility 
DANAE/KLOE-2 planned could start data tak- 
ing in 2010 and provide a further indispensable 
step in the improvement program on a e s(Mz) 

and a^ ad . Note that a 1% measurement in the 

__ 

range 1 to 2 GeV would reduce the error of J7J| to 

A «had(- s o) data = 0.007364 ± 0.000060. 

A challenging long term project for lattice 
QCD is the direct determination of the Adler 



function from first principles. First attempts 
were made in the past |55j . however, the non- 
pcrturbative part is largely dominated by the 
proper inclusion of the tttt — > p resonance. The 
dominant low energy part is extremely sensitive 
to mass and width of the p [7] , and to get them 
correct requires to perform simulations at the 
physical parameters of coupling and light quark 
masses. 

Acknowledgments 

Special thanks go to A. Denig, S. Eidclman, 
W. Kluge, K. Monig, S. Muller, F. Nguyen, G. 
Pancheri and G. Venanzoni for numerous stimu- 
lating discussions. I am very grateful to H. Meyer, 

0. Tarasov and G. Venanzoni for carefully read- 
ing the manuscript. I also thank Frascati Na- 
tional Laboratory and the KLOE group for the 
kind hospitality extended to me. This work was 
supported in part by DFG Sonderforschungsbere- 
ich Transregio 9-03 and by the European Commu- 
nity's Human Potential Program under contract 
HPRN-CT-2002-00311 EURIDICE and the TARI 
Program under contract PJI3-CT-2004-506078. 

REFERENCES 

1. F. Jegerlehncr, in 2nd ECFA/DESY Study 
1998-2001, p. 1851, |hep-ph/0105283| 
hep-ph/0308117; J. Phys. G 29 (2003) 
101. 

2. J. A. Aguilar-Saavedra et al., 
[ECFA/DESY L C Physics Working Group], 
hep-ph/0106315 

3. S. Eidelman, F. Jegerlehner, Z. Phys. C 67 
(1995) 585. 

4. R. Alemany, M. Davier, A. Hocker, Eur. Phys. 
J. C 2 (1998) 123. 

5. M. Davier, S. Eidelman, A. Hocker, Z. Zhang, 
Eur. Phys. J. C 27 (2003) 497; ibid. 31 (2003) 
503 

6. K. Hagiwara, A. D. Martin, D. Nomura, 
T. Teubner, Phys. Lett. B 557 (2003) 69; 
Phys. Rev. D 69 (2004) 093003. 

7. S. Ghozzi, F. Jegerlehner, Phys. Lett. B 583 
(2004) 222. 

8. V. V. Ezhela, S. B. Lugovsky, O. V. Zenin, 



|hep-ph/0312114| 

9. J. F. de Troconiz, F. J. Yndurain, Phys. Rev. 
D 71 (2005) 073008. 

10. H. Burkhardt, B. Pietrzyk, Phys. Rev. D 72 
(2005) 057501. 

11. S. Spagnolo, Eur. Phys. J. C 6 (1999) 
637; A. B. Arbuzov, E. A. Ku- 
raev, N. P. Merenkov, L. Trentadue, 
JHEP 12 (1998) 009; S. Binner, J. H. Kiihn, 
K. Melnikov, Phys. Lett. B 459 (1999) 279; 
M. I. Konchatnij, N. P. Merenkov, JETP 
Lett. 69 (1999) 811; V. A. Khoze et al, Eur. 
Phys. J. C 18 (2001) 481; 

12. M. Benayoun, S. I. Eidelman, 
V. N. Ivanchenko, Z. K. Silagadze, Mod. 
Phys. Lett. A 14 (1999) 2605. 

13. H. Czyz, J. H. Kiihn, Eur. Phys. J. C 18 
(2001) 497; G. Rodrigo, H. Czyz, J. H. Kiihn, 
M. Szopa, Eur. Phys. J. C 24 (2002) 71. 

14. A. B. Arbuzov, D. Haidt, C. Mat- 
teuzzi, M. Paganoni, L. Trentadue, Eur. 
Phys. J. C 34 (2004) 267; L. Trentadue, 
hep-ph/0608076. 

15. A. Sirlin, R. Zucchini, Phys. Rev. Lett. 57 
(1986) 1994; W. J. Marciano, A. Sirlin, Phys. 
Rev. Lett. 56 (1986) 22. 

16. D. H. Wilkinson, J. Phys. G 29 (2003) 189. 

17. Y. M. Shatunov [VEPP-2000 Team Collab- 
oration], Status of the VEPP-2000 collider 
project, eConf C0309101 (2003) WEPL004. 

18. F. Ambrosino et al., Prospects for e + e~ 
physics at Frascati between the $ and the tp, 
hep-ex/0603056 

19. S. G. Gorishnii, A. L. Kataev, S. A. Larin, 
Phys. Lett. B 259 (1991) 144; L. R. Surgu- 
ladze, M. A. Samuel, Phys. Rev. Lett. 66 
(1991) 560 [Erratum-ibid. 66 (1991) 2416]; 
K. G. Chetyrkin, Phys. Lett. B 391 (1997) 
402. 

20. K. G. Chetyrkin, J. H. Kiihn, Phys. Lett. B 
342 (1995) 356; K. G. Chetyrkin, R. V. Har- 
lander, J. H. Kiihn, Nucl. Phys. B 586 (2000) 
56 [Erratum-ibid. B 634 (2002) 413]. 

21. R. V. Harlander, M. Steinhauser, Comput. 
Phys. Commun. 153 (2003) 244. 

22. A. E. Blinov et al. [MD-1 Collaboration], Z. 
Phys. C 70 (1996) 31 

23. J. Z. Bai et al. [BES Collaboration], Phys. 



Rev. Lett. 84 (2000) 594; Phys. Rev. Lett. 
88 (2002) 101802. 

24. R. R. Akhmetshin et al. [CMD-2 Collabora- 
tion], Phys. Lett. B 578 (2004) 285; Phys. 
Lett. B 527 (2002) 161. 

25. A. Aloisio et al. [KLOE Collaboration], Phys. 
Lett. B 606 (2005) 12. 

26. M. N. Achasov et al. [SND Collaboration], 
hep-ex/0605013 

27. V. M. Aulchenko et al. [CMD-2 Collabora- 
tion], JETP Lett. 82 (2005) 743 [Pisma Zh. 
Eksp. Teor. Fiz. 82 (2005) 841]. 

28. M. N. Achasov et al. [SND Collaboration], 
JETP Lett. 68, 573 (1998); JETP Lett. 72, 
282 (2000) [Pisma Zh. Eksp. Teor. Fiz. 72, 
411 (2000)]; Eur. Phys. J. C 12 (2000) 25; 
Phys. Lett. B 486 (2000) 29; AIP Conf. Proc. 
619 (2002) 721; Phys. Lett. B 537 (2002) 201; 
Phys. Lett. B 559 (2003) 171; Phys. Rev. D 
74 (2006) 014016. 

29. R. R. Akhmetshin et al. [CMD-2 Collabora- 
tion], Phys. Lett. B 460, 242 (1999); 494 
(2000) 26; 509 (2001) 217; 562 (2003) 173; 
580 (2004) 119; 605 (2005) 26. 

30. R. R. Akhmetshin et al. [CMD-2 Collabora- 
tion], Phys. Lett. B 489 (2000) 125. 

31. S. I. Eidelman [CMD-2 and SND Collabora- 
tions], Nucl. Phys. Proc. Suppl. 144 (2005) 
223. 

32. M. N. Achasov et al. [SND Collaboration], 
Phys. Rev. D 66 (2002) 032001; ibid. 68 
(2003) 052006. 

33. M. N. Achasov et al. [SND Collaboration], J. 
Exp. Theor. Phys. 96 (2003) 789 [Zh. Eksp. 
Teor. Fiz. 123 (2003) 899]. 

34. R. R. Akhmetshin et al. [CMD-2 Collabora- 
tion], Phys. Lett. B 595 (2004) 101. 

35. R. R. Akhmetshin et al. [CMD-2 Collabora- 
tion], Phys. Lett. B 466 (1999) 385 [Erratum- 
ibid. B 508 (2001) 217]. 

36. M. Ablikim et al. [BES Collaboration], Phys. 
Lett. B 630 (2005) 14. 

37. T. K. Pedlar et al. [CLEO Collaboration], 
Phys. Rev. Lett. 95 (2005) 261803. 

38. M. R. Whalley, J. Phys. G 29 (2003) Al 

39. B. Aubert et al. [BABAR Collaboration], 
Phys. Rev. D 70 (2004) 072004; 71 (2005) 
052001; 73 (2006) 012005; 73 (2006) 052003. 



40. B. A. Shwartz [Belle Collaboration], Nucl. 
Phys. Proc. Suppl. 144 (2005) 245. 

41. A. Hoefer, J. Gluza, F. Jegcrlehner, Eur. 
Phys. J. C 24 (2002) 51; J. Gluza, A. Hoe- 
fer, S. Jadach, F. Jegerlehner, Eur. Phys. J. 
C 28 (2003) 261. 

42. K. K. Seth, Phys. Rev. D 72 (2005) 017501. 

43. A. Osterheld et al. [Crystal Ball Collabora- 
tion], SLAC-PUB-4160 

44. S. Eidelman, F. Jegerlehner, A. L. Kataev, 
O. Vcretin, Phys. Lett. B 454 (1999) 369. 

45. K. A. Milton, I. L. Solovtsov, O. P. Solovtsova, 
Mod. Phys. Lett. A 21 (2006) 1355. 

46. F. Jegerlehner, in "Radiative Corrections", 
ed. J. Sola, World Scientific, Singapore, 1999; 
hep-ph/9901386| 

47. S. Bethke, arXiv:hep-ex/0606035. 

48. J. H. Kiihn, M. Steinhauser, Nucl. Phys. 
B 619 (2001) 588 [Erratum-ibid. B 640 
(2002) 415]; JHEP 0210 (2002) 018; J. Er- 
ler, M. x. Luo, Phys. Lett. B 558 (2003) 
125; M. Eidemuller, Phys. Rev. D 67 (2003) 
113002; G. Corcella, A. H. Hoang, Phys. Lett. 
B 554 (2003) 133. 

49. J. Rolf, S. Sint [ALPHA Collaboration], 
JHEP 0212 (2002) 007; D. Becirevic, V. Lu- 
bicz, G. Martinelli, Phys. Lett. B 524 
(2002) 115; G. M. de Divitiis, M. Guagnelli, 
R. Petronzio, N. Tantalo, F. Palombi, Nucl. 
Phys. B 675 (2003) 309. 

50. W. M. Yao et al. [Particle Data Group Col- 
laboration], J. Phys. G 33 (2006) 1. 

51. V. Gimcncz, L. Giusti, G. Martinelli, F. Ra- 
puano, JHEP 0003 (2000) 018; J. Heitger, 
R. Sommcr [ALPHA collaboration], Nucl. 
Phys. Proc. Suppl. 106 (2002) 358; F. Di 
Renzo, L. Scorzato, JHEP 0411 (2004) 
036; C. McNeile, C. Michael, G. Thompson 
[UKQCD Collaboration], Phys. Lett. B 600 
(2004) 77. 

52. R. Sommer, arXiv:hep-ph/0607088 

53. S. A. Dytman [CLEO Collaboration], Nucl. 
Phys. Proc. Suppl. 13 1 (2004 ) 32. 

54. F. A. Harris, |arX~iv:physics/ 0606059 

55. M. Gockeler et al., Nucl. Phys. Proc. Suppl. 
94 (2001) 571; T. Blum, Phys. Rev. Lett. 
91 (2003) 052001; C. Aubin, T. Blum, 
hcp-lat/0608011. 



